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Abstract 

^ ' We investigate some properties of Bethe-Salpeter wave functions in integrable 

models. In particular we illustrate the application of the operator product expansion 
in determining the short distance behavior. The energy dependence of the potentials 
obtained from such wave functions is studied, and further we discuss the (limited) 
phenomenological significance of zero-energy potentials. 



1 Introduction 



In a recent paper p] which has received general recognition |2j, Ishii, Hatsuda and 
one of the present authors (S. A) have presented results on the nucleon-nucleon (NN) 
potential from first principle lattice computations |3|4j . The results qualitatively 
resemble phenomenological NN potentials which are employed in nuclear physics. 
The force at medium to long range (r > 1.2fm) is attractive; this feature which is 
essential for the existence of bound states of nuclei (e.g. the deuteron) has long 
well been understood in terms of pion and other heavier meson exchange. At short 
distances a characteristic repulsive core is produced [1] by the QCD dynamics, but 
this feature has not yet found a simpler theoretical explanation. 

Intuitively the short distance behavior in QCD is encoded in operator prod- 
uct expansions (OPE). However wave functions and potentials in the framework 
of relativistic quantum field theory are notoriously "flexible" concepts. There are 
infinitely many definitions depending on the interpolating fields chosen and thus 
the universality of the short distance behavior extracted from one particular cho- 
sen wave function is not a priori clear. The wave function discussed in ref. [1] is a 
Bethe-Salpeter (BS) wave function with a particular nucleon interpolating field of 
lowest dimension. The phenomenological success of the results gives rise to the hope 
that one is on the "correct track" , however there remain many theoretical questions 
and refinements in the measurements to be made. For example the results are still 
in the quenched approximation, lattice artifacts must be studied in more detail, the 
dependence of the results on the interpolating field must be examined and the very 
definition of a potential via a BS wave function must be better understood. 

It is our hope that studies of BS wave functions in integrable models in two 
dimensions will give us more insight into such questions. As an aside here we note 
that the methods used in ref. [1] were partially motivated by a method proposed to 
measure phase shifts in a two-dimensional model [5]. In this paper we investigate 
BS wave functions in the Ising model and in the 0(3) non-linear sigma-model in 
two space-time dimensions. 

In a remarkable paper, Fonseca and Zamolodchikov [6j obtained an exact ex- 
pression for the BS wave function of the Ising field theory. In Sect. 2 we study its 
properties; in particular we can see at which distances the short distance behavior 
expected from the OPE sets in. We also point out that a zero energy potential 
defined from the BS wave function has a non-trivial form and may be a concept 
which may have a wider domain of applicability. 

We further examine how the wave function in the Ising model is built from the 
contributions of the intermediate particle states. We find that intermediate states 
involving a relatively low number of particles give a good approximation down to 
quite short distances. This study was performed because in other integrable models 
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the exact wave function is not (yet) known and the only analytical methods available 
are intermediate state approximations, the OPE (renormalized perturbation theory 
at short distances) and 1/n expansions. As an example in Sect. 3 we study the 
(asymptotically free) 0(3) sigma model in two dimensions. 

Various technical details are relegated to appendices and in Sect. 4 we make 
some concluding remarks. 

2 The two— dimensional Ising model in the scaling limit 

In this section we will discuss properties of BS wave functions in the two-dimensional 
Ising field theory, but before introducing these we first briefly describe the theory 
and establish some notations and conventions. 

The theoretical insight which is to be gained from the 2-d Ising model seems 
inexhaustible. In 1976 Wu, McCoy, Tracey and Barouch [7j showed that the model 
(at zero external field) has a continuum limit as one approaches the critical point. 
This relativistic quantum field theory, called the Ising field theory, describes on-shell 
free particles of mass M > 0. We denote the corresponding one-particle states with 
momentum p = M(cosh#,sinh#) by \9), with state normalization 



The continuum limit of the spin field cr(x) is an interpolating field for this 
particle; we chose the normalization 



Although the theory has an alternative representation in terms of free (fermion) 
fields, the spin field is not a free field; nevertheless there is a wealth of information 
on its correlation functions. In ref. [7] it was shown that the two-point function of 
cr(x) satisfies the Painleve III equation. More explicitly defining the vacuum 2-point 
function of a(x) and also that of the corresponding disorder variable /j,(x) at equal 
times (which suffices for our considerations) by 





(0\a(x)\e)=e- 



p = M(cosh 9, sinh 9) . 




G(r) 
G(r) 



(0\a((0, Xl ))a(0)\0) 
<0|/i((0,xi))/i(0)|0) 



(2.3) 
(2.4) 



where r = Mx\. Then for the sum and difference 




(2.5) 



1 which is local with respect to itself but non-local wrt a(x) 
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one has 



G±(r) = e xW/2 e ±¥»(r)/2 t ( 2 _g) 



where the functions Xj V 3 obey the equations 

-W(r)]' = ~sm>(r)), (2.7) 
r 2 

-[rx'W]' = ^[l-ch(2^(r))]. (2.8) 
r 2 

For definiteness we are considering the theory obtained by taking the continuum 
limit from the symmetric phase, where fJ,(x) has a non-vanishing vacuum expectation 
value B 

The short and long distance behaviors of the functions 92, x are summarized in 
Appendix A, and from these it follows that for small r > (for the field normalization 
given in (I2T2D ). 

G(r) ~C x r- 1 / 4 + 0(r 3 / 4 lnr), (2.9) 

(where the constant C x is given in (|A.4j) .) exhibiting the well known anomalous 
dimension of a{x). For large distances r > the correlation function falls exponen- 
tially: 

G(r) ~ -Jl [1 + O^ 1 )] . (2.10) 
Vo7rr 

The results on the 2-point function were subsequently derived in other ways 
(see e.g. [8] and [9]). A very elegant derivation recently presented by Fonseca and 
Zamolodchikov [6] is based on local conservation laws of the doubled Ising field 
theory. 



2.1 Bethe Salpeter wave functions 

In their remarkable paper Fonseca and Zamolodchikov [6 J showed that their methods 
also lead to exact results for a larger class of correlation functions. In particular they 
obtained exact expressions for the BS wave functions for 2-particle in-states 

*(r,0) = ma{{^ Xl ))ame,-e)' m , (2.11) 
*(r,0) = i(Q\^ Xl ))^)\6,-ey a . (2.12) 

In fact Fonseca and Zamolodchikov consider the wave functions for general space- 
time arguments of the fields but here we restrict ourselves to equal times. Without 

2 i.e. the fields a, \i in Fonseca and Zamolodchikov [6] are interchanged with respect to ours. Also 
the field normalization differs. 
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loss of generality we can consider the rapidity 9 > and in the following consider 
only r > since locality (and parity invariance) imply 

¥(r,0) = ¥(-r,0). (2.13) 

For the sum and difference 

*± = $±\P, (2.14) 
Fonseca and Zamolodchikov [6j obtain 

<Mr,0) = [e- e $ ± (r,0) 2 -e 9 $ T (r,0) 2 ] , (2.15) 
where G± are defined in (12. 5h and <I>-|- H satisfy the coupled equations 

*' ± (r, 0) = ^sh(^(r) ± #)c& T (r, 6) , (2.16) 
and the boundary conditions for small r are 

e xM/2$ ± (r, 0) ~ v ^ : C' x e ±e / 2 r 1 / 4 [l + 0(r 2 lnr)] . (2.17) 
In terms of these functions the BS wave function VP is given by 



eXM/2 



$+(r, fl) 2 cosh ( ^ - J - (r, #) 2 cosh ( ^ + , 



(2.18) 



For short distances r it has the expansion 

*(r,0) ~ * as M) + 0(r 7 / 4 ), (2.19) 
^ as (r,#) = 27rC x r 3/4 sinh0, (2.20) 

which is as expected from the known operator expansion of the product of cr-fields 

a((0, xi))<t(0) ~ G(r) + cr 3 / 4 £(0) + .... (2.21) 

where 5 (a;) is the mass operator of dimension 1. 

The coupled differential equations (j2.7jl . (|2.8j) . and <|2. 16j> for $±, <p, x with their 
known boundary conditions at r = can be easily solved numerically. In Fig. Q] we 
depict the wave function \l/(r, 6) for various rapidities, illustrating the early set in 
of the long distance sinusoidal behavior. Fig. [2] shows the wave function divided by 
sinh# so that the leading short distance behaviors are the same, (see (j2.20jl ): once 
this renormalization is done there is rather little remaining variation with the energy 
for r < 0.5; moreover the leading OPE behavior dominates up to r ~ 0.2. 



3 which are the functions \P± in the notation of ref. [6], 
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Figure 1: The Ising BS wave function *(r,0) for 6 = 1.0 (dotted), 6 = 0.6 (dot-dashed), 6 = 0.3 
(dashed), and the zero-energy wave function £(r) defined in Appendix B (solid). 




Figure 2: A renormalized Ising BS wave function ^(r, 0)/sinh(0) for for 6 = (top), 6 — 
0.3,0.6, 1.0 (bottom solid curve). The leading short distance OPE behavior 27rC x r 3/4 (see (J272QJ) 
is given by the dashed curve. 
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Figure 3: The Ising BS potentials (multiplied by r 2 ) r 2 V e (r) for 8 = 1.0 (dotted), 8 = 0.6 
(dot-dashed), 8 = 0.3 (dashed), and 8 = (solid). 



2.2 BS Potentials 

From the BS wave function one can define a rapidity-dependent potential by 

tt''(r,fl)+ S inh 2 flV&(r,fl) 
W •= ^ • (2-22) 

This definition is a direct analogy to that of energy dependent NN potentials made 
in ref. pQ. The hope is that for low energies and for the distances relevant for 
phenomenology the potentials are only mildly energy dependent @. It is such an 
ansatz which seems to qualitatively apply in the NN case. We can investigate this 
question for the Ising field theory and find indeed only moderate variations in a 
reasonable range of parameters, as is illustrated in Fig. [3j The potential for 8 = 1.0 
becomes singular already at r ~ 2.681 where the corresponding wave function has 
its first node. Of course the physics in the Ising model is vastly different from the 
NN case, in particular in the Ising field theory there are no bound states. 

The paper ref. [3] describes some ideas to obtain a local energy independent 
potential from the BS wave functions and this will hopefully be elucidated in our 
next paper [TO]. Here we remark that a natural candidate for a potential of limited 
phenomenological relevance, as we will discuss below, is the zero-energy potential, 



4 Of course this excludes distances near and beyond the point where ^(r, 8) has its first zero and 
hence at which Ve(r) is singular. 
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obtained as the zero energy limit Vo(r) of (|2,22p . For purposes of numerical evalua- 
tion in the Ising model this can be expressed in terms of (p, x (see Appendix A). A 
plot of this potential is included in Fig. [3j The asymptotic behaviors are analytically 
directly obtained from the formulae in Appendix A. From the large r behavior of 
the zero energy wave function 

* (r) := Jim [^(r, 9)] ~ 2r + (l - H + . . . ) , (2.23) 

follows the leading large r asymptotics of the zero energy potential: 

i.e. the potential falls exponentially to zero from above. On the other hand for small 
r using (|2.20p we have 

Vb(r)~-~- ( 2 - 25 ) 
lb r z 

Since this potential is (classically) strongly attractive close to the origin the 
question of possible bound states naturally emerges. This would be fatal for the 
hope that the zero-energy potential is at all relevant for the Ising field theory. In 
Appendix A we show that indeed there are no bound states because (|2.25p is not 
attractive enough in the quantum theory. 

In Appendix B we consider the zero-energy potential in a slightly more general 
context. There we show that it reproduces the correct scattering length, which 
parameterizes the leading low momentum behavior of the phase shift. However in 
general it does not yield the exact next-to-leading behavior (although it may in some 
cases be a good approximation to it). 

2.3 Intermediate particle state approximations to ^>(r, 9) 

In the Ising field theory we are, as discussed above, fortunate to have exact partial 
differential equations for the BS wave functions. However for most other integrable 
models this is not (yet) the case, and we have to resort to approximations in order 
to obtain quantitative results. One approach is to compute contributions from inter- 
mediate states from knowledge of the form factors. For the two-point function this 
approximation has been investigated in ref. |llj . and there the contributions of only 
a few states is found to approximate the exact result down to very small distances 
where the OPE can be applied. 

In the Hilbert space of in-states defined by the spin field the S-matrix operator 
is given by 

S = (_i)N(N-i)/2 } (2 .26) 
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where N is the particle number operator. An energy independent phase is not 
observable in a scattering experiment; however the non-trivial S-matrix reflects the 
fact that cr(x) is not a free field. 

Given knowledge of the S-matrix and assuming general properties such as an- 
alyticity and crossing symmetry (together with some additional technical assump- 
tions) it was argued in ref. [12] that generalized form factors of the spin field are 
given by 



out / 



h,...,0Mo)\e t+1 ,...,e n r 

= (2 ,)(n-i)/2 Yl r(|^-%|) J] T( V _g) II T(\0k-0i\), 

l<*<j<* \<r<t<s<n ^ r s ' t<k<l<n 

(2.27) 

with n an odd positive integer. V denotes the principle part and 

T(z) = tanh|. (2.28) 
Sandwiching a complete set of states 

' - ioxoi + § /_: 1: . . . ah*,.. . . , m <*»> 

(where s stands for in or out) between the fields, \P can be expressed as a sum over 
s-particle contributions 

oo 

^(r,e)=J2^2 S -i(r,9). (2.30) 

s=l 

Starting with the 1-particle contribution we have 

^i(r,9) = -^-T(20)pi(r,e). (2.31) 



with 



Now use 



to obtain 



p lM )=/ i e ie ^ w -^- Try (2.32, 



OO 



V 1 

2vri<5(0) + — — , (2.33) 



Pl=P? '+P? + P?, (2.34) 
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where the superscripts denote the number of delta functions. So 

/•oo 

pf\r,d) = (2vri) 2 / d9 1 e irsMi 5(e 1 -e)5(e 1 + e) 

J — oo 



-2n 2 5(9) . 



(2.35) 
(2.36) 



Next 



2m 



2-ki 



d0i e 



irshO 



S(8i - 6) i <5(gi + 0) 
T{9 l +0 + ie) T(9 1 -9 + ie 



c 



irshO 



-irsh9 



T(20 + ie) T{20-ie) 
sin(rsh#) 



An 2 5{9) - 4vr 



T(20) 



Finally 



Pi (r,< 



A() pirshBi 

oo 1 T(0i - + ie)T{0 x + + ie) 



(2.37) 
(2.38) 

(2.39) 



(2.40) 



Shifting the contour to the line parallel to the real axis with imaginary part in/ 2 
(and observing that the contribution from the contours parallel to the imaginary 
axis at infinity is zero for we get 



(0)/ 

Pi [r, 



Summarizing we have 

*i(r,0) 



dze~ r ' 



rchz 



ch9 + ishz 
ch9 — ishz 



2K (r) -4ch 2 6» / dz 



-rchz 



2ir 



ch 2 8 + sh 2 z 



4vrsin(rsh6l) - T(20)pf\r, 8) 



(2.41) 
(2.42) 

(2.43) 



The plane wave part is as expected for a two particle S-matrix equal to —1; and 
Pi\r,0) decays exponentially as r — > oo. However v I / i(r, 9) diverges logarithmically 
as r — > 0: 



*i(r,0)~-T(20) pnr + /(0) + O(r)] , 

7T 



(2.44) 



which is very different from the short distance behavior of the exact wave function 
given in (j2T20|) . 

The contribution from the 3-particle states can be computed similarly. Here we 
just note that all contributions vanish exponentially as r — > oo (some only as e _r due 
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to disconnected contributions). We have numerically computed these contributions 
for various rapidities and as a typical result we give the results for rapidity 6 = 0.3 
in Table Q] where we compare the approximations to the exact wave function. We 
observe that whereas (for this rapidity) the 1-particle approximation fails quite badly 
at r = 0.1, addition of the 3-particle contribution already makes the agreement much 
better at this distance and already here the asymptotic formula (|2.20p which can 
be derived from the OPE sets in. Addition of the 5-particle intermediate states 
would of course improve the agreement to smaller distances as illustrated in the 
0(3) o"-model in the next section. 



r 


*i(r,0.3) 


tf 3 (r,0.3) 


[*i + *3](r, 


0.3) 


*(r,0.3) 


tf M (r,0.3) 


10.0 


1.92482e - 


1 


1.39916e - 


6 


1.92484e - 


- 1 


1.92484e - 


1 




5.0 


1.99793 




2.52957e - 


4 


1.99818 




1.99818 






4.0 


1.87759 




7.18688e - 


4 


1.87831 




1.87831 






3.0 


1.58539 




2.03727e - 


3 


1.58743 




1.58743 






2.0 


1.14970 




5.73985e - 


3 


1.15544 




1.15544 






1.0 


6.12374e - 


1 


1.63733e - 


2 


6.28747e - 


- 1 


6.28748e - 


1 




0.1 


2.60053e - 


4 


9.40535e - 


2 


9.43136e - 


- 2 


9.45618e - 


2 


9.227e - 2 


0.01 


-2.43697e - 


1 


2.50800e - 


1 


7.10370e- 


-3 


1.64495e - 


2 


1.641e - 2 


0.001 


-4.59999e - 


1 


3.85735e - 


1 


-7.42642e - 


- 2 


2.91863e - 


3 


2.918e - 3 


0.0001 














5.18897e - 


4 


5.189e - 4 



Table 1: ^(r, 6) the exact wave function, and ^i(r, 9), ^(r, 6) the 1- and 3-particle contributions, 
and the leading short distance behavior ty as (r8) given in (|2.20|l . for 6 = 0.3. 
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3 BS wave functions of the 0(3) a— model 



In this section we will give a quantitative discussion of BS wave functions and their 
associated potentials in the 0(3) non-linear sigma model in two dimensions. The 
0(n) sigma model has long served as a favorite laboratory for testing ideas concern- 
ing asymptotically free theories |13|14|15|16|17j . Unfortunately there is in this case 
no exact expression known for correlation functions of any local operators. However 
for the case n = 3 the multi-particle form factors (FF) (defined in Eq. (|3.17p ) can 
be obtained recursively (although they become extremely complicated for higher 
particle states), and so one can often obtain excellent approximations to correlation 
functions in a wide range of energies by saturation with a low number of intermediate 
states. 

The spectrum is considered to contain an 0(n) vector multiplet of particles of 
mass M and to have no bound states. The two-particle S-matrix established by 
Zamolodchikov and Zamolodchikov [18] is given by: 

2 

S ab ;cd(P) = E Si(f3)Pi(ab\cd) , (3.1) 

1=0 

where (3 is the rapidity difference of the incoming particles and Pi are "isospin" 
projectors given byd 

P (ab\cd) = -5 ab 5 cd , (3.2) 
n 

1 1 

Pi(ab\cd) = -5 ac 5 b d - ^ ad S bc , (3.3) 

111 

P 2 (ab\cd) = -5 ac 5 bd + -5 a d5bc &ab&cd , (3.4) 

2 2 n 



and 

S x {p) = -(-l) 7 e 2 ^ (/3) . (3.5) 
For the special case n = 3 the phase shifts are simply given by 

(5b 09) = arctan^ , (3.7) 
W) = (3-8) 



5 Note J2 e ,f Pi{ab\ef)Pj{ef\cd) = 5/ jP^atycd) , and for n = 3 one has £ s P z (ab\ab) =21 + 1. 
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We define BS wave functions as in the last section by 

Vab;cd(xi,0) = (0|a a (0,xiy fe (0,0)|c,M,-#> in , 9>0. (3.9) 

The spin field a a (x) is an interpolating field for the massive particle and we fix the 
normalization by 

(0|<r a (0)M) = 5 ab . (3.10) 
Translation invariance and locality implies 

Vab;cd(-Xl,e) = ^ba;dc{xi,0). (3.11) 
Introducing isospin components Tj for all tensors T^.^ : 

Ti = ^-jPjiablcd^cd , (3.12) 

Eq. (|3.1ip now implies 

* I (-x 1 ,9) = (-l) I * I (x 1 ,9). (3.13) 

As before the BS wave function can be expressed as an expansion over s-particle 
contributions: 

*«6;cd(*l,0)= E **,cdM)> ( 3 - 14 ) 
s odd 

and these can further be organized in contributions having specific large distance 
behavior. Since the computation is rather technical we relegate the details to Ap- 
pendix C and here just summarize the result. Firstly it is convenient to introduce 
the modified wave function 

* /(r ' 9) = t^h8 e ~ i5im * l{r/M > 0) ' (3 ' 15) 

since as shown in Appendix C ^ becomes real (for real arguments). This WF has a 
(large distance) expansion which is naturally expressed in the form 

h(r,6)= £ {A^\r,9) + B ( [ m \r,9)} , (3.16) 

m odd 

where A { ™\r, 9) ~ 0(e- mr ) and s} m) (r, 9) ~ 0(e - ( m - 1 ) r ) for large r. Their explicit 
expressions are given in Eqs. (|C.3ip . (|C.32|) and involve integrals over products of 
the generalized form factors 

(0|<7 a (0)|6i,/3i; . . . ;b s ,f3 s )' m = ^..^(fr, ■ ■ ■ ,&) , Pi > fa > • • • (3.17) 
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-II ' ' ' 1 1 1 

0.5 1 1.5 

r 

Figure 4: Contributions to the 0(3) wave function in the 7 = channel for 9 = 0.3. The dotted 
curve is B£\r,Q.S), the dashed curve is A.^ ^(z - , 0.3), the dot-dashed curve is B^ 3) (r,0.3), and the 
solid curve is the sum of the first 5 contributions in the long distance expansion. 

From the connectivity properties of the matrix elements built of these form factors 
it follows that the s-particle contribution (s > 3) contributes not only to and 
but also to A( s ~ 2 \ The leading term is 

2 

B { j L \r,6) = — sin {rsinh0 + 5/(20)} . (3.18) 

tanh 

We have numerically computed the first 5 terms in the expansion for the three 
isospin values (in the case of 0(3)) and a range of small (< 1.0) rapidities. As for 
the case of the Ising model, we find that inclusion of sufficient number of terms in 
the long distance expansion gives a good description of the full wave function down 
to quite small distances. This is illustrated in Tables [3HS in Appendix D where we 
give results for 9 = 0.3. Although the individual terms diverge at short distances 
it seems that their sums are tending to zero in each channel see Figs SHSl Exactly 
how the limit is reached can however not be read off from this approximation and 
we require a detailed OPE analysis which we will present in the next subsection. 

As in the case of the Ising models the BS wave functions and their corresponding 
potentials are only weakly varying with the rapidity (for moderate rapidities) in the 
short to intermediate distance range. 

Note for the zero energy WF we have 

B ( I 1) (r,0)=2(-l) I {r + a I } , (3.19) 
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where 

1 1 2 . 

a = , a\ = — , a2 = — . (3.20) 

7T 7T 7T 

For large r the potentials fall exponentially to zero; in particular for the zero energy 
potentials: 



1 9-TT 3 

= ^—[l-^r-' + Oir-^e-r, (3.21) 



1 / ?7T 3 

V x {r) = -_y-^[l-0.056r- 1 + O(r- 2 )]e- r , (3.22) 

1 /27T 3 

F 2 (r) = [1- 1.570^ + 0(r- 2 )]e- r . (3.23) 

3.1 Operator product expansion 

Consider the operator product 

Df{y) = P^ablcd) (* c (y)a d (0)) , (3.24) 

for general n > 3. Using asymptotic freedom, we can show that its leading short 
distance expansion is of the form 

D?(y)^a(\y\)t ab (0) + ..., 

D?(y)*(3(\y\)y»jf(0) + ..., 

sab , ~ > (3.25) 

D${y) « — { 7 o(|y|) + |y| 2 7i(|y|)e(o) + ^72(1^1)^(0)) + . . . 



Here t is a traceless iso-tensor operator of dimension 0, is the Noether current 
and and Tn V is the (Lorentz) trace and traceless part of the energy-momentum 
tensor T^ u : 

e = t^, v = % v + ^e, (3.26) 

r/oo = —T/ii = 1. Finally the leading short distance behavior of the functions 
a(|y|), 7j(|y|) appearing in (|3.25j) can be computed in the framework of renor- 
malized perturbation theory as discussed below. 

Sandwiching the operator product between the vacuum and a two-particle state 
we have 

(0\Df(y)\c, 9; d, -9) in = Pi(ab\cd)^ j(r, 6) , (3.27) 
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where r = M|y|. 

Let us recall |19| the two-particle form factors of the operators occurring in the 
above short distance expansion: 



(O\t ab (O)\c,a;d,0) in = (a - - m) tanh ^^y^J ffe(a&|cd), ( 3 - 28 ) 
(O\jf(O)\c,a;d,0) [n = -i^e^tpia - 0)P 1 (ab\cd) , (3.29) 



(O\T, u (O)\a,a;b,0r = - q,q u ) ^ _ P) . (3.30) 

Here 

q° = M(cosh a + cosh 0), q 1 = M(sinha + sinh/5) , (3.31) 



eoi = -eio = 1 and 



Eq. (|3.28p fixes the normalization of t ab , which is otherwise undetermined. 

Using these form factors, the leading short distance expansion of the isospin 
invariant wave functions are given as 

^ 2 (r,0) ps a(r)(20-z7r)tanh0, 
, , „. „, itt 2 sinh 6 tanh 6 „ 

W r/?(r) 26(i«-0) {20 ~ m) ' (3.33) 

*oM)^(^) 2 [2 7l (r)-72(r)]- ^ ° 



49 (iw — 9) 

which in terms of the redefined (real) field in (|3.15p read: 



~ , a . a . .^ 2 sinh0 /vr 2 + A6 2 

*!(r, (?) « r 0(r)—^—yJ -^-^ , (3.34) 



2 ... „ sinh 26* 1 



# Q (r,0) ~ (r7r) z [2 7l (r)- 72 (r)] 



Vvr 2 + 2 
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We now outline the information which can be gained from perturbative field 
theory; for any undefined notation we refer the reader to [20]. We start with the 
short distance expansion 

Af (y) = \Pi(ab\cd)S c (y)S d (0) ^ifj% 05 y)S? )o6 (0) + ..., (3.35) 

where j3^ ab are (bare) local operators and the are coefficient functions (which 
can, in principle, be calculated in perturbation theory). The operator product (|3.24j) 
differs by a (non-perturbative) rescaling from the renormalized version of (|3.35j) : 



n-'Df(y) = Af (R] (y) = P J (o6|cd)5^ ) (y)5f fl) (0) 

,(w)o 
h(R) 



£a>Wv)s$?(o), (3 - 36) 



written in terms of renormalized operators S^j^ b and finite coefficient functions kj . 
The latter satisfies the renormalization group (RG) equation 

[v + i{g)+ if\g)} 4 W) =0, (3.37) 



where the RG 7-function 



lf\g)=1^9 2 + ... (3.38) 



is related to the operator renormalization constant (in dimensional regularization) 
by 

Z M = : _ 7^ + _ _ _ ^ {3 3g) 

corresponding to the operator B^> ab . 

In particular, for I = 2 we have only one operator 

B ab = f ab = }_( gagb _ ± § ab\ (3 4 q) 

% V n J 

which has renormalization constant [20] 

t 

ITS 

and coefficient function 



Z 2 = Z f = ! + — + ... (3.41) 



K 2 (g ,y) = l + O(g 2 ). (3.42) 
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Similarly, for / = 1 we have 



Bt = Jt, (3-43) 



with Z\ = 1 and 



^f(5o,2/) = -^ + 0( 5o 2 ). (3.44) 



Finally for / = we have the two operators 

4^ = 76- B$? = ^ V (3-45) 

with = 1 (o» = 1, 2) and coefficient functions 

4 1] = C w \y\ 2 + O( 5o 2 ), K^ v = -lyV + O( 5o 2 ) . (3.46) 

Since vanishes at tree level, the corresponding numerical coefficient C\q can only 
be determined by a one-loop calculation. 

The equation (|3.37|) can be solved by standard RG methods. Introducing the 
running coupling function A(r) as the solution of 

— - + xhiA(r) = -In r, (3.47) 
A(r) 

with x = V( n — 2) we find 

a(r) » a A x {l + O(A)} , 

/3(r)«-^ n (27r X Ar x {l + 0(A)} , 

1 (3.48) 
71 (r) « C w D n (27r X \r x {1 + 0(A)} , 

72 (r)«-^„(27r X A)-^{l + 0(A)} . 

The constant occurring in the coefficient a cannot be calculated since we do not 
know the relative normalization of the operator with respect to t ab whose nor- 
malization is fixed non-perturbatively by (|3.28[) . We also do not know the numerical 
value of the non-perturbative constant D n for general n. However, we do know [20] 
D3 = 4/-7T and this enables us to write for n = 3 (and, for simplicity, at zero energy): 

^o(r,0)«^(c 10 + £){l + 4 1) A + ...}, 

*i(r, 0) ^{l + c^A + ...}, ( 3 - 49 ) 
§ 2 (r,0) «a ^A{l + 4 1) A + ...} . 
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Here the 0(A) (and higher) corrections can in principle be calculated in higher 
orders of perturbation theory. The number Cio can also be obtained by a one-loop 
calculation. However, as explained above, the coefficient «o cannot be calculated by 
presently available methods. 

Fortunately, the overall normalization cancels from the potential defined by 



V7W 



yi(r,0) 
*/(r,0) 



(3.50) 



and we find 



V (r) ~ 


r 2 { 


V x (r) ~ 




r 2 i 


V 2 (r) ~ 


r 2 1 



1-±A + 0(A 2 ) 



1 + 



1 



l-c 



(i) 



l-c. 



(1) 



A + 0(A 2 )} , 
A + 0(A 2 )} . 



(3.51) 



Note that for I = we can calculate the first correction in A(r) without knowledge 
of the 0(A) correction in f|3.49[) . 



r 


0.01 


0.05 


0.1 


0.2 


A 


0.154 


0.222 


0.280 


0.393 



Table 2: Values of A(r). 

In Fig. [7] we plot r 2 times the potentials in the I = 1,2 channels obtained 
from the sum of the first 5 leading terms in the long distance (LD) expansions, 
together with the leading behavior (|3,5ip obtained from the OPE. They are plotted 
with respect to the variable A defined in (|3.47p (with x = l)j i n Table [2] we give 
some pairs of values (r, A(r)). We also plot the curve —A + 2A 2 to illustrate that 
("quite reasonable") higher order expressions in the OPE expansion could be found 
to make smooth meetings with the LD approximations. We think that the 5-term LD 
approximation is accurate down to A ~ 0.2 (which already corresponds to quite short 
distances) in the 1 = 1 channel and even to smaller distances in the 1 = 2 channel. 
This can be monitored by studying the stability of successive LD approximations, 
including only 3 terms, 4 terms and 5 terms respectively. This is illustrated in Fig. [8] 
for 1 = 1. Alternativey one can appreciate the situation by inspecting Table [6] in 
Appendix D where we give the double derivatives of the separate contributions times 
r 2 in the various channels. 

Fig. [9] shows the zero-energy potential in the 1 = channel. Here the LD breaks 
down already larger distances, in fact there is no stability in the sense described 
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0.05 



- 



-0.05 ■ 



-0.1 



-0.15 



-0.2 i 



0.1 



0.2 0.3 
lambda 



0.4 



Figure 7: "Long distance approximation" to r 2 Vi(r) (solid) and r 2 V2(r) (dashed). The lower 
dotted line is the leading short distance behavior —A and the upper dotted line is —A + 2A 2 . 



-0.1 



-%1 



0.15 



0.2 0.25 0.3 

lambda 



0.35 



0.4 



Figure 8: Successive LD approximations to r 2 Vi(r); The dotted, dashed and solid lines correspond 
to approximations using 3,4,5 terms repectively. 
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1 - 



0.5 - 




lambda 



Figure 9: LD approximation to r 2 Vo(r) (solid). The upper dotted line is the leading short distance 
behavior and the lower dotted line is 2 — 3A — 8A 2 . 

above even at r = 0.1. The figure however suggests that our 5-term approximation 
may still be quite good there, but it would need computation of higher order terms 
to confirm this. Never the less it is plausible that the approximation joins smoothly 
to the OPE expansion where again higher order terms are required to improve the 
quantitative picture. 
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4 Conclusions 



In this paper we have investigated BS wave functions in integrable models; the Ising 
and 0(3) cr-models. We have seen that potentials derived from them are rather 
slowly varying with energy in the short and intermediate distance range M\x\ <~ 1. 
We have also discussed the relevance of the zero-energy potential and its phenomeno- 
logical limitations. In these models we have found that a good approximation to the 
wave functions can be obtained by combining a long distance expansion (from con- 
tributions of intermediate states) and a short distance expansion from the OPE. It 
would be instructive to study some other examples in particular models with bound 
states. 

Given a BS wave function constructed from a particular choice of local fields, it 
is clear that its short distance behavior can be obtained by an analysis of the OPE 
expansion. Further it follows from naive dimensional analysis that the BS potentials 
derived from (most) wave functions that vanish at the origin will behave as \x\~ 2 
(modified by logarithms). The overall sign in 3-dimensions indicates its attractive 
or repulsive nature. How this sign may depend on the particular interpolating field 
remains an important question. 

In a sequel paper [10] we plan to include a more general discussion on potentials 
obtained from BS wave functions. We will also present OPE predictions for the short 
distance behaviors of BS wave functions (and the resulting short distance behavior 
of the potentials) for the pion-pion and nucleon-nucleon cases in QCD, for some 
choice of the interpolating fields.. 
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Appendix A 

A.l Asymptotic behaviors of the functions ip, x 

Note here we always take r > 0. The short distance behaviors of 93, \ are given by 

e-^W ~ ~\rVL{r) [l + 0(r 4 )] , (A.l) 
e x(r) „ -2C^Q(r) [1 + 0(r 2 )] , (A.2) 



where 



0(r) = ln(fcr) , A; = -e 7 , (A.3) 



where 7 is the Euler-Mascheroni constant (= 0.57721566490 . . . ), and 

C x = 2- 7/6 A" 3 e 1/4 = 0.27119012339 . . . , (A.4) 
where A is Glaisher's constant: 

A = exp|^-C'(-l)| = 1.282427... (A.5) 

Next for long distances [6]: 

<p{r) = -K (r) + O(e- 3r ), (A.6) 

7T 

X (r) = -21n2-^ [r{Ko(r) 2 -^i(r) 2 } + i^ (r)^i(r)] +0(e- 4r ). (A.7) 
Recall the modified Bessel function 

POO 

K (r) = I dze~ rchz , (A.8) 

behaves for small r as 







^o(r) = -lny- 7 + 0(r 2 ), (A.J 



and for large r, 



Ao(» 



1-^ + ... 



(A.10) 



and Ki(r) = -K' (r) . 
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A. 2 Zero energy potential 

Expanding 3>±(r, 9) for small 6: 

$ ± (r,0) = $ o (r)±0$i(r) + . 

and denning 

/i(r) = $oW 2 , 

v{r) = 2$„(r)$ 1 (r), 

we have 

fjf(r) = fj,(r) sinh <^(r) , 

i/(r) = |u(r) cosh <^(r) . 



Then it follows for small 9: 
tf(r,0) ~ 26>e x(r)/2 



z/(r) cosh ^ - — /x(r) sinh ^ ^ 



+ ... 



and 



Vo(r) := lim 



¥"(r,0) 



X '(r)(-+ X '(r))-^(r) 



+ 



F(r) 



(A.ll) 



(A.12) 
(A.13) 



(A.14) 
(A.15) 



(A.16) 



v{r) cosh — /x(r) sinh 



.(A.17) 



with 



Fir) = \ {xWto + <p"(r)} |i/(r) sinh ^ - M ( r )cosh^} 



+/x(r) {x'(r) + sinhy?(r)} cosh 



y>(r) 



(A.18) 



A. 3 Absence of bound states 

A bound state would have negative energy — E and would be a solution of the 
Schrodinger equation 



-<£»" + V <S> = -E$ . 



(A.19) 



with large r asymptotics 



(A.20) 
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i.e. would be normalizable. 

It is important to impose the right boundary conditions at the origin. It turns 
out [21|22j that the correct boundary condition is 

$(r)~r 3/4 . (A.21) 

In other words, the other, more singular solution of the Schrodinger equation, which 
would behave like 

$(r) ~ r 1 / 4 (A.22) 

is not allowed (the Hamilton operator would not be self-adjoint). 
Using this information it is easy to derive the formula 

rr 

tf' (r)$(r) - V (r)&(r) = ~ E I dx <S>(x)^o(x) . (A.23) 

J o 

The left hand side vanishes for large r and we know that the zero energy wave func- 
tion ^>o(r) (|2,23p is positive everywhere. Therefore the bound state wave function 

has to change sign somewhere. We denote this point by ro and assume <&(r) is 
positive between the origin and this point. Then we have 

rro 

-*o(ro)$'(ro) = ~E / dx $(x)V (x) , (A.24) 
J o 

implying that 

$'(r )>0, (A.25) 
but this is obviously a contradiction, so $(r) cannot exist. 
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Appendix B. Scattering length and effective range 

Let us consider the BS wave function 

$(r,fc) = (0|a(0,r)a(0,0)|2) (B.l) 

in a 2-dimensional model, where k is the wave number of the 2-particle state. Its 
large r asymptotics 



$(r,Jfe) ~ <p{k)sm{kr + <5(A;)} (B.2) 



can be used [5] to read off the physical phase shift 5{k). This has a low energy 
expansion of the form 

6(k) = -ak + fk 3 + --- (B.3) 

It is convenient to introduce a new normalization for the wave function so that for 
large r we have 

sin < kr + 5(k) > 

$(r, k) ^ >- . (B.4) 

k 

The zero energy wave function 

l(r) = $(r,0), (B.5) 

has long distance asymptotics 

£(r)~ r -& + ■■■ (B.6) 

It can be used to define the zero energy potential 

i"(r) 

We can now study the solution of the Schrodinger equation with this potential: 
-<(r, k) + U(r)if)(r, k) = k 2 ip(r, k) , (B.8) 
where the solution is fixed by requiring the following asymptotic behavior: 

„ , , sin {kr + 5(k)} ,„ , 

tf>(r,k) 1 ' ( R9 ) 

where for low energy 

8(k) = -ak + fk 3 + • • • (B.10) 
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Since the zero energy wave function is the same as before, 

^(r,0)=£(r), (B.ll) 

we have 

a = a, (B.12) 

i.e. the scattering length corresponding to the zero energy potential is equal to the 
physical one. 

Using the Schrodinger equation (|B.8P we can write 

A {V'(r, k 2 Mr, fa) - V'(r, fa)^{r, k 2 )} = (kf - k 2 2 )^(r, fa^r, k 2 ) , (B.13) 
which can be integrated to give 

i)\RM)^{RM)-i>\RM)^{RM) = {kl-k 2 2 ) / driP(r,h)?Jj(r, k 2 ) ■ (B.14) 

Jo 

Now we define S(r, k) by 

S{r,k) ;= ^{fa' + ^)} > (B .15) 



which satisfies 



{kl-k 2 2 ) I drS(r,k l )S{r,k 2 ) 
Jo 

= S'(R, k 2 )S(R, fa) - S'(R, fa)S(R, k 2 ) + S'(0, fa)S(0, k 2 ) - S'(0, k 2 )S(0, fa) . 

(B.16) 

Taking the difference between (|B.16P and (|B.14p in the R — > oo limit gives 



(k\ - k\) I dr {S{r, fa)S(r, k 2 ) - ^(r, fa)^(r, k 2 )} 







sm5(k 2 ) sm5(ki] 

COsd{fa) COS0(K2j 



ko fa 



Now we take k 2 — > and fcx = fc in the above formula and find 

/>oo 



fc 2 / dr {(r - a)5(r, fc) - ^(r, fe)i(r)} = -a cos <5(Jfe) - sin ^ fc ) . (B.17) 



/o 
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Expanding this exact formula we get 

,3 



f = Y~ B ' (B.18) 



where 



/•oo 

B= I dr {{r-af - i 2 (r)} . (B.19) 
Jo 

In scattering theory one usually introduces the effective range: 

P = ^ 2 , (B.20) 

which, together with the scattering length a, gives a two-parameter description of 
low energy scattering. 
In the Ising model 

a = a = 0, f = 0, (B.21) 

and p cannot be defined. On the other handd, 

/■OO 

f = -B= dr {£ 2 (r) - r 2 } ~ 0.263. (B.22) 
Jo 

We see that the zero energy potential does (in general) not reproduce the low energy 
expansion of the true phase shift beyond leading order. 



6 The fact that / > follows directly from the property that £(r) > r for all r. 
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Appendix C. s— particle contributions to the 0(3) wave functions 



The matrix elements built from the form factors (|3.17p have a particular connectivity 
structure (see e.g. ref. |23j). from which one infers that the s-particle contribution 
can be written as a sum of three terms with D delta-functions involving the rapidity 
variables, D = 0, 1, 2: 



*S^i.fl) = E*S2 ) («i.f). (ci) 



D=0 

We start with contributions without delta functions: 



* ( 6 ( d( a; i! 9 ) = I d/3i . . . d/3 s ^ l[sinhft+ .,, +sinh& 



//9i>...>A ( 47r ) s (C.2) 

X FbL.tsiPll ■ ■ ■ ■>Ps)J : b s ...b 1 cd(Ps, ■■■,$1,0, -0) , 

where 

& = Pi + iTr-ie. (C.3) 

We first note that here the integrand is a totally symmetric function of the (real) in- 
tegration variables Pi. Thus we can extend the integration from the original domain 
Pi >■■■> Ps to H s . After that we can shift the integration contour by defining 

Pi = cxi — ivr/2, cti real . (C4) 

This latter step works for negative x\ only, and for this reason from now on we take 

Mxi = -r, r>0, (C.5) 

and set M = 1 in the following. We will use (|3.13p later to get the wave function 
for positive x\. We get 

^■!S(-r,0) = ^ cd (r,e), (C.6) 

where 

W <„ a\ — 1 f°° • • • ^ a £ p -r[coshQ!i+-+cosha s ] 

x ^6i...6 8 ( a l>- ' ' ^ a s)^b s ...b 1 cd{° l s,- ■ . ,ax,e - —,-e - —) . 



s! (4vr) , c ^ 



This behaves as 0(e _sr ) for large r. 

Next we discuss contributions with two delta-functions. Let us discuss first the 
case s = 3. Here we see using the FF axioms that the three terms have the same 
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analytic form, the only difference is that the range of the integration variable j3 is 
different, namely 



-6 > p, 6 > (3 > -9, (3 > 9 (C.8) 

for the three terms. This means that the sum of the three contributions can be 
simply written as the same integral with the integration extending from — oo to oo. 
Similar considerations work for general s and we get 

V( l ' )_ l>...> ft (4vr)(- 2 ) (C.9) 

x J r cdb- i ...b a {9i ~^,^3, • • • iPs)^b s ...b 3 0s, ■■■)/%)• 

We can again extend the integrations to IR/ S_2 ) and then shift the integration 
contours: 



H,(a){ 2 )(_ r Q\ _ 1 f°° dQ 3 • ' • da s -r-[cosha 3 +-+eosha s 

^ ' j ( S -2)!7_ 00 (4tt)(-2) 



Z7T zvr 

x ^6 3 ...6 S ( + -J ' _6> + y ' a3 ' • • • ' ^^...bsi^s, ■ ■ ■ » «3) 



(C.IO) 



Finally we use the relation (expressing parity invariance) 

. . . , e s ) = rt.. bl (-e s , . . . , , (en) 

and see that ([C.lOp can be expressed with (|C.7|) : 

^■S(-^) = F^(r,e). (C.12) 

The last group of integrals is with one delta-function. In the s = 3 case we can 
group the six contributions into two groups of three. The first one contains integrals 
proportional to e lXl sinh e , whereas the integrals in the second one are proportional to 
g-j^isinhe anc j a ^ go con t a i n the factor S c d-..(^)- Again, the three terms in the first 
group are of the same analytic form and correspond to the domains 

6>f3>(3', 0>9>0', p>(3'>9, (C.13) 

for the two integration variables /3,/3'. The sum of the three integrals is simply an 
integral over j3 > /3' . Generalizing to arbitrary s the first group gives 

ab ' cd ( ' j 4>...>ft (4*)<-u" ^ (C.14) 

x Fcb 2 ...b s (0'P 2 > ■ ■ ■ >Ps)Fb s ...b 2 d0s, ■■■,02, -6) ■ 
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Performing the usual operations we get 

0) = e- irS ' mh9 9 ( :U r > °) + ^ rS ' mhd Sc d -M29)g^, dl (r, -9) , (C.15) 

where 



„(») <V m _ 1 f°° d«2 • • • d«s ^- r [cosha2+-+coshq a 



(s-iyj-oo (4^)( s - 1 ) 

ITT 

x Kb 2 ...bS e + y, «2, • • • , ««)^6 J ,...6 a d( Q! «> • • • > "2, -6 1 - y) 



(C.16) 



It is easy to see that this behaves as 0(e~( s_1 ) r ) for large r. 

Adding up all the contributions and using (|3.13p we can write for the isospin 
components: 

^M) = 2(-i) 7 F / (s) M) 

s odd (C.17) 
+ £ {{-lYe- irsinhe g\ s) (r^O) + e irsinhd S I (29)g i I s) (r, -9)^ , 

s odd 

where 5/ are the isospin invariant S-matrix amplitudes in (|3.ip . Equation (|C.17p 
is a large distance expansion and we now study the first few terms (up to the ones 
behaving 0(e _3r ) for large r). 
For 0(1) we get 

gf\r,9) = l. (C.18) 

For 0(e _r ) we have 

F 7 (1) (r,0) = ^(29) f°° dae- rcosh Xa + y-0)V(a + y + 9) Pl (a,9), (C.19) 



where 

Po (a,9) = -49 -2m, p x (a, 9) = m - 2a , pz(a,9) = 29 - 2m. (C.20) 
Here we have used the representation (for s > 2): 

KM^---^s) = 7r 3 ^/ 2 g a bl ^ bs (9 1 ,...,9 s ) H m-Oj), (C21) 

l<i<j<s 

with i/j(9) defined in (13,321) and where g b b (9\, . . . , 9 S ) is a polynomial function of 
the rapidities. Some explicit expressions of the latter are given in |17j . 
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Next we have 

-4 f-oo 



4 /"OO 

9?\r,9) = da 2 da 3 e- r ( cosha2+cosha3 )/ l/ (a 2 ,«3^) 

x V(«2 - 03)^(03 - a 2 )ip{6 + y - a 2 )ip(9 + y - a 3 ) (C.22) 

x $(6 + y + a 2 )iP(e + — + a 3 ) , 

where /i/ is the quadratic polynomial 

hi(a 2) a s ,e) = 2I + l Pi(ab\cd)g% 2bi {6 + —,a 2 ,az)gl 3b2d (a ?> ,a2,-0 - —) . 

(C.23) 

Finally 

if } M) = |^(20) y da 1 da 2 dQ 3 e-'-( coshQl + cosh ^+ cosha3 )u; / (« 1) a 2 ,a3^) 
3 

x ^(a fe + y + + y ~ ^) II ~ - «fc) 

fc=i ~ " fc</ 



(C.24) 



with 



u;/(ai,a 2 ,a 3 ,<9) = 2J 1 Pj(o&|cd)gg l&263 (ai, a 2 , a 3 ) 

X 9b 3 b 2 b lC d( a 3, «2, ai, 6» - y , -0 - y , 



(C.25) 



Let us now study the phase of the WF. We start from the relation 

{J^... 6s (0i, ...,8 S )Y= K...b s (-^ (C26) 

expressing the fact that the FF is a real analytic function (which is a consequence 
of CPT symmetry, but can also be proven directly). Using this in (|C.7|) we find the 
relation 

= S cd , yx {-2e)F^ xy {r,6) , (C.27) 

which gives 

{F 7 (s) (r,60}* = (-l) / 5 / (-20)F / (s) (r,0) = -e-^^if^r, 6) . (C.28) 
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For (|C.16P we simply get 

{gth<r>°)Y = 9%*(r-0), (C.29) 

and 

{s^M)}* = gf\r,-9). (C.30) 

Using (fCTfj) and the relations ([028|) . ([O30]) we see that $ defined in (f37T5]) is real. 
Further the functions occurring in its long distance expansion f)3. 16H are given by 

^ m)M) = ^ke { - 1)Ie ~ lSli29)F i m){r ' e) > (c - 31) 

and 

B^\r,9) = -^(-lYluiL-^^e-^^g^ir^)} , (C.32) 
tanh u I J 

with F( m ) , </ m ) given in Eqs. (|C.7|) . and (|C.16P respectively. 
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Appendix D. 0(3) a model tables 



r 


^ lj (r,0.3) 


< } (r,0.3) 


^ 3) (r,0.3) 




^ 5, (r,0.3) 


sum 


10.0 


1.30735 


1.06995e - 5 


-8.766e - 12 


8.0801e - 21 


-3.1e- 28 


1.30736 


5.0 


6.79501 


1.96768e - 3 


6.5182e - 7 


2.6482e - 13 


2.98e - 18 


6.79698 


4.0 


6.18820 


5.6340e - 3 


8.8594e - 6 


1.0298e - 11 


1.40e - 15 


6.19384 


3.0 


5.01196 


1.61265e - 2 


1.1513e-4 


4.5395e - 10 


5.51e- 13 


5.02820 


2.0 


3.37453 


4.55362e - 2 


1.5472e - 3 


2.4111e - 8 


2.86e - 10 


3.42161 


1.0 


1.42658 


0.116843 


2.4630e - 2 


1.6800e - 6 


3.34e - 7 


1.56806 


0.1 


-0.444244 


-0.162087 


0.666238 


-1.3290e - 3 


5.355e - 3 


0.063932 


0.01 


-0.631821 


-0.930030 


1.48854 


-4.1388e - 2 


0.11430 


-0.000397 


0.001 


-0.650555 


-1.53211 


1.93324 


-0.21387 


0.44076 


-0.02252 


0.0001 


-0.652428 


-1.93408 


2.12862 


-0.54207 


0.9129 


-0.0870 



Table 3: s = 1, 3 and _B (5) contributions to 0(3) isospin wave functions for 6 = 0.3. 
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T 


a l \r, U.6) 


A^\r, 0.3) 


u^hr n "i\ 
n l [r, u.6) 


Aj 3J (r,0.3) 




U.OJ 


sum 


10.0 


-0.00197770 


5.8085e - 


6 


— 1.84e - 12 


-7.4550e - 


21 


-l.Oe 


- 29 


-0.00197711 


5.0 


-6.85843 


1.1598e - 


3 


5.9610e - 8 


-2.7141e - 


13 


1.41e - 


- 19 


-6.85727 


4.0 


-6.63613 


3.4459e - 


3 


1.0439e - 6 


-1.1041e - 


11 


5.16e - 


- 17 


-6.63268 


3.0 


—5 80318 


1.0444e - 


2 


1 5554e — 5 


-5.209e - 10 


1.82e - 


- 14 


—5 79272 


2.0 


-4.43624 


3.2680e - 


2 


2.2998e - 4 


-3.1147e - 


8 


8.44e - 


- 12 


-4.40333 


1.0 


-2.66108 


0.108668 




3.9490e - 3 


-2.9371e - 


6 


8.31e - 


- 9 


-2.54846 


0.1 


-0.848874 


0.368994 




0.102435 


-2.7406e - 


4 


7.407e 


- 5 


-0.37764 


0.01 


-0.661862 


0.445499 




0.157496 


6.0312e - 3 




1.528e 


- 3 


-0.05131 


0.001 


-0.643131 


0.524720 




7.8661e - 2 


2.3673e - 2 




9.049e 


- 3 


-0.00703 


0.0001 


-0.641257 


0.621496 




-0.054472 


0.046377 




0.024426 


-0.00343 



Table 4: As in Table [3] but for isospin 1. 



r 


4>,0.3) 




^ 3j (r,0.3) 


Af{r, 0.3) 


^ 5, (r,0.3) 


sum 


10.0 


-0.632946 


1.0557e - 5 


-2.03e - 12 


1.211e- 20 


-1.9e - 30 


-0.632935 


5.0 


6.79781 


1.9415e - 3 


5.5835e - 8 


4.6314e - 13 


-1.06e- 19 


6.79975 


4.0 


6.77335 


5.5592e - 3 


1.1078e - 6 


1.9250e - 11 


-1.35e- 17 


6.77891 


3.0 


6.12563 


1.5912e - 2 


1.7834e - 5 


9.3827e - 10 


-7.34e - 16 


6.14156 


2.0 


4.91424 


4.4931e - 2 


2.8628e - 4 


5.9387e - 8 


1.19e - 12 


4.95945 


1.0 


3.25064 


0.115291 


5.6063e - 3 


6.4031e - 6 


3.37e - 9 


3.37155 


0.1 


1.49266 


-0.159933 


0.235546 


2.5843e - 3 


9.673e - 5 


1.57095 


0.01 


1.30846 


-0.917675 


0.680703 


1.10920e - 2 


2.688e - 4 


1.08285 


0.001 


1.28998 


-1.51175 


1.03886 


3.0193e - 2 


-0.010478 


0.83681 


0.0001 


1.28814 


-1.90839 


1.27147 


7.1185e - 2 


-0.037702 


0.68470 



Table 5: As in Table E] but for isospm Z. 
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r 


r 2 A ( > (r,0) 


r 2 5^ j (r,0) 


o , 3)" 


(r,0) 




(r,0) 


r 2 V (r) 


5.0 


5.312e-2 


2.60e - 4 


8.86e - 


11 


2.26e - 


14 


5.70e - 3 


2.0 


0.16415 


5.273e - 2 


1.66e - 


6 


7.99e - 


8 


6.36e - 2 


1.0 


-2.875e - 2 


0.24430 


3.10e- 


5 


2.71e - 


5 


0.1436 


0.2 


-0.39396 


0.47923 


— 2.88e 


— 3 


1.035e - 


- 2 


0.54* 


0.1 


-0.39508 


0.41926 


-9.72e 


- 3 


3.086e - 


- 2 


** 


r 


r 2 ^ j (r,0) 


r 2 Bf> (r,0) 




(r,0) 


r 2 Bf> 


(r,0) 


r 2 Vi(r) 


5.0 


3.439e - 2 


4.66e - 5 


-8.98e 


- 11 


5.94e - 


16 


-3.24e - 3 


2.0 


0.18153 


9.74e - 3 


-2.27e 


-6 


1.77e - 


9 


-4.15e - 2 


1.0 


0.17885 


4.503e - 2 


-7.28e 


- 5 


5.03e - 


7 


-8.87e - 2 


n 1 

U.l 


i.io4e — z 


4.1UUe — Z 


1 K A ^ 

i.o4e — 


Q 
O 


o.zoe — 


A 

4 


n i /ion 
— U. 14yU 


0.05 


1.275e - 2 


1.135e - 2 


3.761e 


-3 


7.17e - 


4 


-0.141* 


0.01 


3.256e - 2 


-3.763e - 2 


7.74e - 


3 


2.92e - 


3 


** 


r 


r^'V, 0) 


r 2 B^ r '(r, 0) 


r 2 Af 


(r,0) 


r 2 Bf 


(r,0) 


r 2 y 2 (r) 


5.0 


5.312e-2 


5.52e - 5 


1.55e- 


10 


3.16e- 


16 


4.72e - 3 


2.0 


0.16415 


1.35e - 2 


4.51e- 


6 


1.46e - 


9 


3.34e - 2 


1.0 


-2.875e - 2 


7.439e - 2 


1.77e - 


4 


5.93e - 


7 


1.35e - 2 


0.1 


-0.39508 


0.23011 


3.99e - 


3 


6.39e - 


4 


-0.1039 


0.01 


-0.27014 


0.16771 


6.21e- 


3 


-3.22e 


- 3 


-9.3e- 2* 



Table 6: Double derivatives of the contributions to 0(3) wave functions for 6 = in all isospin 
channels. A double star ** indicates that there is no stability and a single star * indicates a ~ 10% 
variation between successive approximations. 
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